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Accuracy of the mean activity coefficient expression (Hansen-Vieillefosse-Belloni equation), valid within the 
hypernetted chain (HNC) approximation, was tested in a wide concentration range against new Monte Carlo 
(MC) data for +l:-l and +2:-2 primitive model electrolytes. The expression has an advantage that the 
excess chemical potential can be obtained directly without invoking the time consuming Gibbs-Duhem cal- 
culation. We found the HNC results for the mean activity coefficient to be in good agreement with the ma- 
chine calculations performed for the same model. In addition, the thermodynamic consistency of the HNC 
approximation was tested. The mean activity coefficients, calculated via the Gibbs-Duhem equation, seem to 
follow the MC data slightly better than the Hansen-Vieillefosse-Belloni expression. For completeness of the 
calculation, the HNC excess internal energies and osmotic coefficients are also presented. These results are 
compared with the calculations based on other theories commonly used to describe electrolyte solutions, such 
as the mean spherical approximation, Pitzer's extension of the Debye-HCickel theory, and the Debye-Hiickel 
limiting law. 
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1. Introduction 

Knowledge of the physico-chemical properties of electrolyte solutions contributes toward better 
understanding of the processes in biology, chemical sciences, and various technologies. It is of no 
surprise that electrolyte solutions are, after a hundred years, still extensively studied in search of 
theories providing quantitative description of these systems. The advances in modeling, together 
with the review of experimental data for electrolyte solutions, are summarized in books and the- 
matic papers (see, for example, references ITjiSj). Although it has been established that only the 
models that explicitly include solvent are capable of providing a proper microscopical description 
of electrolyte solutions jT] , one can describe thermodynamic properties of these systems on the 
McMillan-Mayer level of approximation by a proper choice of ion size parameters [5J [TO] • In this 
paper we will restrict ourselves to the latter case. 

Seminal efforts to understand the properties of electrolyte solutions were made by Debye and 
Hiickel (DH) [11] . The DH theory is based on a linearized version of the Poisson-Boltzmann equa- 
tion which yields important insights into dilute electrolyte solutions. Unfortunately, this approach 
can only be used to quantitatively describe very dilute solutions of size symmetric -j-l:— 1 elec- 
trolytes [T5] . The deficiencies of the Debye and Hiickel theory were analyzed in many contributions 
(see, e.g., references pi 151 [TOMl6| l. Various continuations were suggested to extend the range of 
applicability of the theory, among these the one proposed by Pitzer [17] is very useful for en- 
gineering purposes. This approach is simple, fully analytical, and capable of semi-quantitatively 

* Dedicated to Professor Yura V. Kalyuzhnyi on the occasion of his 60"" birthday. 
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describing the osmotic coefficient behavior of +1: — 1 electrolyte up to the concentration close to 
1.0 moldm-3 [T7]. 

Theoretical extension, correcting for the deficiency of the Debye and Hiickel theory, is provided 
by the modified Poisson-Boltzmann approach |14 | ll8lf2T| . For example, the individual activity co- 
efficients of pure electrolytes were calculated by Molero et al. (cf. reference [22 ) and for ternary 
systems containing neutral hard spheres recently by Outhwaite et al. [23 . The theory yields excel- 
lent agreement with computer simulations. 

Another approach to calculate the properties of solutions involve a class of integral equation the- 
ories based on the Ornstein-Zernike (OZ) integral equation |24|E^ and approximate closures. The 
so-called mean spherical approximation (MSA) was solved analytically for the primitive model elec- 
trolytes, with the thermodynamic quantities written in a closed form [26 -33 . Through the energy 
route, one obtains the osmotic and activity coefficients that are in good agreement with Monte Carlo 
computer simulations for -fl: — 1 electrolytes [12; . Due to its analytical nature, the MSA theory 
provides useful insights and is accordingly extensively used (see, e.g., references P^l - HO] !. Another 
approximate but very robust closure of the Ornstein-Zernike theory was proposed by Kovalenko 
and Hirata (KH) fUlfi^]. By changing the radii of the ions j32lH5H15] and/or including the associa- 
tion between unlike ions |46l 147] these methods provide good fits to experimental data. An integral 
equation approximation that is widely used in describing the thermodynamics of symmetric, as 
well as highly asymmetric electrolytes, is the hypernetted-chain (HNC) closure [T21 1211 HHl HI] • 

An advantage of the OZ hypernetted-chain approximation is that it yields accurate structural 
description of electrolyte solutions in terms of the pair-distribution functions, gij{r). Once this in- 
formation is known, the standard statistical-mechanical equations, connecting the pair-distribution 
functions with thermodynamic properties, can be applied | ,24j . The reduced excess internal energy 
of the system is readily obtained via the expression 

^^ = fEE— /"..Wft.Mdr, {^,J^+,-) (1.1) 
t j P J 

where j3 = (fee being the Boltzmann constant and T the absolute temperature), pi (pj) are 

the number densities of the species i (j), p is the total number density of the system, Uij{r) the 
interaction pair potential (see sectional), and dr = Anr'^dr. 

The osmotic coefficient, $hnc , can be conveniently calculated using the virial route [24 . For the 
primitive electrolyte, the pair potential contains a discontinuity at the distance of closest approach 
of two ions (see section [2]) . The problem can be surmounted by the introduction of the so-called 
background correlation (referred to as the cavity distribution) function, yij{r) = exp[/3uij{r)]gij{r), 
and the formal division of the Uij{r) to the Coulombic, u^(r), and the hard sphere part. In this 
way, the equation gets the form applicable to the primitive model electrolyte |24| 

W^l-^EE^/^^^^^.Wdr+f ^5:^4,^,(a,), (z,,=+,-). (1.2) 

i j i j 

Here aij denotes the distance of closest approach of two spherical particles i and j, and ufj{r) cx r^^ 
for all r, is the usual Coulomb potential. 

The mean activity coefficient, 7± , can be obtained by integration of the Gibbs-Duhem equa- 
tion [U 

c 

ln7±,GD = *HNC - 1 + y"($HNC - l)dlnc', (1.3) 

c=0 

where c denotes the molar concentration of the electrolyte solution. 

Though the numerical procedure looks straightforward, it may be time consuming: (i) the 
osmotic coefficients need to be obtained for the series of concentrations, and (ii) for electrolyte 



concentrations lower than c = 0.001 moldm , the integral in equation (1.3l may diverge due to the 



numerical problems. This particular problem can be resolved by using the analytical continuation. 
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for example DHLL+B2, at very low electrolyte concentrations fSlP. The mean activity coefficient 
calculation via the Gibbs-Duhem equation becomes particularly cumbersome and time consuming 
in ternary mixtures containing additional solutes |51j . The expressions written above are general 
and can be in principle, apart from the HNC theory, applied to any other approximate closure. 

An alternative way of calculating the activity coefficients, but valid only within the HNC 
approximation, has been proposed by Verlet and Levesque [55]. In its current form the formula was 
written by Hansen and Vieillefosse [53' , and was successfully applied by Belloni [54' for asymmetric 
electrolytes. The expression, here referred to as the Hansen- Vieillefosse-Belloni (HVB) equation, 
reads 

ln7..HVB = -^P,c|f (0) + ^5]p,y {h,{T)[h,,{r)-c,,{r)]}Av, (1.4) 

3 j 

Here 7^ is the individual activity coefficient of the species i (+ or — ), h{r) denotes the total, and 
c(r) the direct correlation function, while c(^)(0) denotes the Fourier transform of the short-range 
part of the direct correlation function at = 0. Further, pj is the number density of species j . 
To our best knowledge, the accuracy of the activity coefficient expression given by HVB equa- 



tion ( 1.4 1 has not been thoroughly tested for the primitive model electrolytes. In a few cases, where 
the Monte Carlo data for the same system were available [55], the expression was found to be in 
good agreement with the "exact" machine calculations. Therefore, the purpose of this work is to 



systematically examine the validity of HVB equation (1.4 1 for the size symmetric and asymmetric 
-|-1:— 1, and size symmetric +2:— 2 electrolytes. Notice that the formula given by equation (1.4) is 
only valid within the HNC approximation and therefore is not generally applicable. In other words, 
different expressions for the mean activity coefficients can be derived consistently with the closure 
conditions used [56'. 

Even so, the HVB equation is very useful because it avoids long and cumbersome evaluation 
of the excess chemical potential of the solute via the Gibbs-Duhem route (see, e.g., reference |51|). 
Notice that the mean activity coefficient of bulk electrolyte is often an input information for 
studying heterogeneous systems, electrical double-layer, or Donnan equilibrium, i.e. wherever the 
bulk electrolyte is in equilibrium with charged surfaces. The results obtained by means of HNC 
approximation, coupled with the HVB equation, are in the present study compared with the new 
Monte Carlo simulations, and with the results of some other electrolyte theories used in describing 
electrolyte solutions. In addition, the internal consistency of the HNC theory is tested with respect 
to the two different routes of the mean activity coefficient calculation. 

The paper is organized as follows: following a brief introduction, the model and methods are 
presented, followed by the results and discussion section. Numerical results are presented in form 
of figures and tables. Conclusions are given at the end. Appendices summarize all the relevant 
equations used in these calculations: Debye-Hiickel theory (appendix A), Pitzer's approach (ap- 
pendix B), and the mean spherical approximation study (appendix C). 



2. Model and methods 

The model most frequently used in describing the electrolyte solutions is called the primitive 
model. In this description, an ion is represented as a hard sphere carrying positive or negative 
charge in its center, while the solution as a whole is treated as a dielectric continuum with a 
relative permittivity Er (McMillan-Mayer level of approximation) of pure solvent at pressure (p) and 
temperature (T) of observation [21] . Despite its simplicity, the model is capable of explaining many 
experimentally determinable properties of real electrolytes [Hi [57], as well as their mixtures [10 . 

The interaction pair potential between two ions of valence Zi {zj), separated by a distance r, 
contains the hard sphere and Coulomb part, u^(r) 

f 00, r < Qij , 

Puij{r) = I (2.1) 
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where = -^{ai + Oj), Oi (o-j) is the diameter of species i (j), and Ab is the Bjerrum length 

Here /3 = l/k^T, where is the Boltzmann constant, T is the absolute temperature, Eq is the 
permittivity in vacuum and eg is the elementary charge. For aqueous solutions at 25 °C studied here, 
the Bjerrum length assumes the value Ab = 7.14 A. Three different models of aqueous electrolyte 
solutions at 25 °C were considered: (a) z+ = — z_ = 1, a+ = a_ = 4.25 A; (b) z+ = — z_ — 1, 
a+ = 5.43 A, a_ = 3.62 A; and (c) z+ = — z_ = 2, a+ = a_ = 4.25 A. These values are often used 
as representatives for simple ions in aqueous solutions. 

2.1. Monte Carlo simulation 

The Monte Carlo calculations were performed in the canonical and grand canonical (GCMC) 
ensemble, using the standard Metropolis sampling algorithm ^58j. Periodic boundary conditions 
with the minimum image (MI) convention, and Ewald summation (ES) method were used in the 
canonical ensemble simulations to minimize the finite sample effects, while GCMC was performed 
only in combination with the Ewald summation. The calculations were carried out with equal 
number of cations and anions, the total number of particles N being between 200 and 1000. After 
an equilibration run of (1 — 10) • 10^ configurations, each production run consisted of (1 — 10) • 10^ 
attempted configurations. From four to ten independent runs were performed for simulation in 
canonical ensemble for each concentration studied. The results given in tables and figures are the 
average values over these runs, after the equilibration runs were discarded. 

While the excess internal energy and osmotic coefficient were obtained as simple canonical 
ensemble averages |5S], the Widom's test particle insertion method |59ll5U] was used in the canonical 
ensemble simulation to calculate the excess chemical potential of the electrolyte in solution 

/3/i± - ln7± = ln(cxp[-/3[/pan.]). (2.2) 

In the expression above, (• • • ) denotes the canonical ensemble average, while C/pair is the interaction 
energy between the system and a non-perturbing test pair of oppositely charged ions, inserted at 
random locations in the system. During the canonical ensemble simulation, for every N (number of 
particles in the simulation box) configuration, a hundred of (Widom's) insertions were attempted, 
to obtain the canonical ensemble average requested by equation ( |2.2[ ). 

The results for the mean chemical potential of the electrolyte solution as obtained by Widom's 
method may, at higher concentrations, depend on the number of particles in the system [61J. For 
this reason, in many cases, especially for more concentrated electrolytes, the canonical ensemble 
simulations were supplemented by the grand canonical ensemble Monte Carlo method. In the latter 
case, the excess chemical potential is obtained directly, without invoking the insertion method |621 

2.2. Hypernetted chain approximation 

The hypernetted chain (HNC) approximation is based on the Ornstein-Zernike (OZ) equation. 
For multi-component mixtures, the OZ equation reads | |24j 

h^,{r) = c,^{r) + Y,Pk I Kk{\r-v'\)ck,{v'W, (2.3) 

k 

where h{r) and c(r) are the total and direct correlation functions, respectively, and the integral 
is of the convolution type. A general closure relation between h{r) and c(r) for the OZ equation 

is HE] 

In [h,j{r) + 1] = -Puijir) + hij{r) - Cij(r) + B,j{r). (2.4) 
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Bij (r) is in literature known as the "bridge graph" and cannot be written as a closed form function 
of the distribution functions h{r) and c(r). Bij{r) are all the graphs in the representation of the 
background correlation function, which are neither series nor parallel graphs. At least two field 
points are bridged, i.e. connected by the Mayer bond. The HNC approximation assumes that these 
graphs mutually cancel and sets Bij{r) to zero. Due to the long-range nature of the Coulomb 
interaction, the set of equations ( |2.3[ ) and (2.4 1 has to be re-normalized before it can be solved 
numerically [IS', . 

The re-normalized form of the integral equation was solved by direct iteration using the fast 
Fourier transform routine on a linear grid with 2^® division points separated by the distance of 
Ar=0.005 A. It should be noted that the results strongly depend on the number of points and 
separation interval Ar. The decisive criteria in our case was the smallest zeroth (electroneutrality 
condition) and second moment condition that should both in theory be equal to zero |65l I66j . 

Equations (1.1 1 and (1.2 1 were used to calculate the excess internal energy and the osmotic 
coefficients, respectively, while the equation (1.4) was used to calculate the activity coefficient. The 
mean activity coefficient, 7-1- , for charge symmetric electrolytes studied here, is finally calculated 
from the expression = 7+7- . In order to check the thermodynamic consistency of the HNC 
approximation, the mean activity coefficients were (for a limited number of cases) calculated with 
the help of the Gibbs-Duhem equation (1.3 1, using the HNC osmotic coefficient (equation (1.2|) 
results. 



3. Results and discussion 

Thermodynamic properties, e.g. excess internal energies, osmotic coefficients, and activity co- 
efficients, are for three examples: (a) z+ = — z_ = 1, a+ = a_ = 4.25 A; (b) = — z_ = 1, 
a+ = 5.43 A, a_ — 3.62 A; and (c) z+ = — z_ — 2, a+ — a- — 4.25 A, as a function of the 
electrolyte molar concentration, presented in tables [T]-[6] and shown in figures [T|]3] The molar con- 
centration of charge symmetric electrolytes, c = c+ = c_ (mol dm^'^) is related to the number 
densities in the following way: pi= Ci Na , where A'a is the Avogadro's number. Basic equations 
for the mean spherical approximation, Pitzer's approach, and the Debye-Hiickel theory, which are 
not in focus of this paper, are given in appendices. 

(a) 2+ — — z_ = 1, a+ — — 4.25 A. 

The excess internal energy, f3E°^/N, chemical potential. In 7-1- , and osmotic coefficient, <&, at 
different concentrations of size symmetric -|-1:— 1 electrolyte are presented in tablesfl] and [2] While 
all the GCMC results were obtained using the Ewald summation technique (table llF, the compar- 
ison between the canonical Monte Carlo results obtained using the periodic boundary conditions 
with minimum image convention and results obtained using Ewald summation are given in table |2] 
For simulations in the canonical ensemble, the Widom method was used to calculate the mean 
activity coefficients. As seen from table[2] for +1: — 1 electrolyte in the concentration range studied 
here (from 0.0001 moldm"^ to 1.5 moldm"^), both ways of accounting for the finite sample effects 
(minimum image, Ewald summation), yield the results, which agree within 5% or (most often) 
better. This appears to be within the sum of numerical uncertainties of separate calculations. 

Figure [1] shows the results presented in tables [l] and |2] (In 7-1- and 1 — i.e. deviations from 
ideality), as well as the results obtained by the frequently used analytical theories (DH approxi- 
mation, Pitzer's approach, MSA theory - see equations in the appendices). In7± (panel (a)) and 
1 — $ (panel (b)) are shown for different methods as a function of square root of the concentration. 
It is clear from this figure that the DH theory quantitatively describes the osmotic and activity 
coefficients of the model electrolyte only up to approximately c = 0.05 moldm"^. Pitzer's continu- 
ation extends the range of the DH approach up to 0.75 moldm"-^. By contrast, the MSA and HNC 
approximations are equally good in the whole concentration range, only at concentrations higher 
than 1 M, the HNC performs slightly better. 



The main purpose of this work is to test the validity of the HVB equation ( 1 .4 1 . First we checked 



the thermodynamic consistency of the HNC approach by comparing the mean activity coefficients 



obtained by the Gibbs-Duhem equation with those obtained directly (equation (1.4|). The former 
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0.4 0.8 

Figure 1. In7± (panel (a)) and 1 — $ (panel (b)) as a function of c^''^ for a +1: — 1 primitive model 
electrolyte, a+ — a- = 4.25 A, Ab ~ 7.14 A. Circles are grand canonical Monte Carlo data, 
crosses represent the canonical ensemble Monte Carlo results (Ewald summation). Theoretical 
values are given by lines: continuous line - HNC, dashed line ~ MSA, dash-dotted line - Pitzer, 
and dotted line - DH. 



results in table |2] are denoted as 7±_gd • One can see that the two sets of results agree very well. 
It is our impression, however, that the mean activity coefficients calculated via the Gibbs-Duhem 
equation are slightly closer to the MC data than those calculated by the HVB equation. The latter 
results for 7-1- seem to be systematically too high, but the differences are small, most likely within 
the experimental uncertainties. 

(b) z+ = = 1, a+ 5.43 A, a_ = 3.62 A. 

The results for size asymmetric +1:— 1 electrolyte are presented in tables [3] and |4] and figure |2] 
Notice that in the DH theory calculations the parameter a was taken to be a = ^(a+ + a_) since 
the DH theory cannot take the size asymmetry of the ions into account. The conclusions are the 
same as for the case (a): the size asymmetry considered here does not dramatically effect the 
performance of these theories in the applied concentration range. 



0.4 




Figure 2. Same as in figure [T] but for a+ = 5.43 A and a_ = 3.62 A. 



(c) z+ = -z_ = 2, a+ = a_ = 4.25 A. 

The excess internal energy, f3E'^^/N, excess chemical potential, ln7± , and osmotic coefficient, 
for different concentrations of size symmetric +2:— 2 electrolyte are presented in tables |5] and [6] 
and figure [3] (only In 7-1- and Interestingly, different treatments of the boundary conditions do 
not effect the results for the excess chemical potential obtained by Widom's method. This is not 
entirely true for the excess internal energy and osmotic coefficient calculations at concentrations 
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above 1 nioldm^'^. The discrepancies may become quite large at concentrations sa 1.5 moldm^'^, 
and the Ewald summation method should be used to obtain correct results. For this reason the 
results of computer simulations, using the minimum image convention at 1.5 moldm^'^, are not 
included in table |6] At such a high concentration of the +2:— 2 electrolyte, the minimum image 
method yields unphysical pair distribution functions, and consequently should not be used for 
accurate simulations. 



-0.6 ■ 



-1.2 ■ 



-2.4 ■ 





Figurs 3. The same as in figure[T] but for the +2:— 2 model electrolyte. 



As seen from figure [3| the differences among different theories become, as expected, more 
pronounced in the case of +2:— 2 electrolyte than for +1:— 1 electrolyte examined before. The 
results presented in this paper suggest, in close agreement with many previous calculations, that 
the HNC approximation accurately describes the thermodynamic properties of the primitive model 
electrolytes in a wide concentration range. 



4. Conclusions 



The expression for the mean activity coefficient valid within the HNC approximation, Hansen- 
Vieillefosse-Belloni equation, was for primitive model +1: — 1 and +2:— 2 electrolytes tested against 
the newly obtained Monte Carlo simulation data. For the sake of completeness, the thermodynamic 
properties of +1: — 1 and +2:— 2 electrolytes, calculated by means of some other theories, often used 
in describing the electrolyte solutions, are also presented. 

Special attention is paid to the numerical accuracy of simulations and the HNC calculations. 
Although the HVB equation (expression (1.4)) has been used before, to our best knowledge, a 
systematic test of its accuracy has not been performed so far. For all the examples studied here, the 
HVB equation, compared to the canonical or grand canonical ensemble Monte Carlo simulations, 
yields accurate results for the mean activity coefficients and is, especially at higher concentration 
and +2:— 2 electrolytes, superior to other approximations examined here. At the same time, it 
appears to be simple in usage and hence represents an excellent tool in describing the excess 
chemical potential of bulk electrolyte solutions. Notice, that this information is needed |55l 170] 
whenever the membrane equilibria involving electrolyte solutions are studied. 
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A. Debye-Hiickel theory 

In the framework of the Debye-Hiickel theory, the expression for the reduced excess internal 
energy reads 

PE^H _ Ab|z+Z-| k 
N ~ 2 l + na' 

where and z_ are valencies of cations and anions, respectively, is the so-called Debye 
screening length defined as = 47rAB ^ Pizf (Ab being the Bjerrum length and pi the number 
density of the ionic species i), and a is the distance of closest approach of two ions (assumed to be 
same for all pairs). 

The osmotic coefficient is given by 

$DH = 1 - ' ' • Kfl , A2 



where a{Ka) is a function defined as |17l IS7] 

<j{x)^ ^(l + x- 21n(l + x)V (A3) 

x-^ \ I + X J 

The mean activity coefficient is calculated via 

2.303A|z+z_|77 ,,,, 

where / is the ionic strength of the solution (/ — O.B^^Cizf , Ci being the molar concentra- 
tion of species i), and A, B are constants containing the absolute temperature and dielectric 
constant of the solvent. For aqueous solutions at 25 °C, A = 0.511 dm'^/^ mol~^/^, and B = 
0.329 • 10-8 cm~i dm3/2 mol-^/^. 

At extreme dilution, the term Ba\/1 in the denominator of equation (A4| becomes negligible 



compared to unity and the equation ( |A4[ ) yields the known Debye-Hiickel's limiting law (DHLL) 

ln7±,DHLL = -2.303A|z+z_|V7. (A5) 
The expression for the osmotic coefficient reduces in this case to |68| 

<I'DHLL = l-^^N+^-|V7. (A6) 



B. Pitzer's approach 

The osmotic coefficient is obtained via the virial route, which yields |17| 



1 



Ab|z+z-| 



6 



1 



-.pa 



Ab|z+z-| 
12a 



1 



(Bl) 



In this expression p ■= p+ + p- \& the total number density of ionic species (for significance of other 
quantities see appendix A) . An important improvement above the Debye-Hiickel theory is the last 
term in equation (Bl ). The term provides the correction to the Debye-Hiickel theory at intermediate 
concentrations and it causes the osmotic coefficient to increase at higher concentrations, exactly 
as observed experimentally. 

The equation for the mean activity coefficient is |17| 



ln7±,P = 



>^b\z+z^\ 
6 



1 



2 Q IAr , , 
-Tipa^ - — L_^ln(l 

3 6a 



$P - 1, 



(B2) 



where $p is given by the equation (Bl I. 
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C. Mean spherical approximation 

Within the MSA closure, one calculates the excess internal energy as [3U] 



f}E\ 



MSA 



a 



(Cl) 



where = 47rAB , A = 1 — ^ piof , and 



n ^ 1 



2A ^ 1 + Toi ' 

i 

p _ 1 piO^iZi 

" ~ f7 ^ r 



1 + ra, 



(C2) 
(C3) 

(C4) 



The electrostatic contribution to the osmotic coefficient is given by 



A 



(C5) 



where p — pi is the total density of the system. The electrostatic contribution to the activity 
coefficient of the species i has the form |5U1 [55] 



in 7,; — — ZiMi — 

47r 4A 



12A 



(C6) 



where 



' ~ 2r(l + rai) 



(C7) 
(C8) 



The hard sphere contribution to the osmotic and activity coefficients follow from the equation 
of state of a mixture of hard spheres. Here we use the Mansoori-Carnahan-Starling-Leland equation 
of state [59]. In this case, the osmotic coefficient reads [M] 



A3 



(C9) 



whereas the hard sphere contribution to the activity coefficient of the species i is 



ln7h^= (^^.-1-^2/3 I - In A 



A2 

A3 



3(1 - ai) + pi + -^{cn - Pi- pi-l) 



9 /3 
^Vz - i^yi - '2' + r] [ -yi - 3?/2 - 4y3 + 1 



22/3 + 4 



(CIO) 
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Notations in equations ( |C9[ ) and ( |C10 ) denote the following 

yi 



and as usual 



j>i 



V2 



^/a^a7 

VkPk x/OiOj 



»J E 



2/3 



77P Ofc 



2/3' 



P V a-iCLj 



p a. 



A, 



Pi ^ -/(HOk 

f]i ^/aja^ 



ik ■ 



]>k 

3 — 2/3 
?7 



6 



rj ai 

2/3 



P 



6 



E^' 



The osmotic and activity coefficients of the primitive model electrolyte solution are then ob- 
tained by summing up the electrostatic and hard sphere contributions 



^MSA 

ln7i,MSA = ln7f + In 7, 



1 , l„-,hs 
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ripi/iMiTMBHi MOfle/ii e/ieKTpo/iiriB. flopiBHJiHHJi peay/ibTariB 
rinep/iaHi4K)roBoro Ha5/ii/i>KeHH5i p,f\9\ Koect3ii4ieHTy aKTi/iBHOcri 
3 flaHMMM MoHTe-Kap/io 

E. ryTieppe3-Ba.n.naflapec!^, M. 71yKiui4^, B. Mi/i.naH-Ma.ncP, B. rpi6ap-JlP, B. B/iaxiPl 

L4eHTp npuK/iaflHoT 4>i3MKM Ta nepeflOBi/ix xexHO/iorii/i, HaLiioHaiibHwti aBXOHOMHUw yHiBepcwTex MexiKO, 
A.C. 1-1010, 76000 Kepexapo, MeKCHKa 

VHiBepcHxex /1ki6/1!Rhw, 4>aKy/ibxex xiwiT xa xiMlHHUx xexHO/ioriw, AcKepneBa 5, SI— 1000, 
JlKi6ji!RHa, Cnosemff 

ToHHicxb Bwpasy fl/iiq KoecjjiLiieHxa cepeflHboT aKxi/iBHOcxi (piBHaHHiq FaHceHa-BiyiJiectDocca-BejiJioHi), cnpaBe- 
fl/iMBoro y rinep/iaHLiKDroBOMy (T/l) Ha6jii/i>KeHHi, nepeBipji/iacn y Lui/ipoKOMy iHxepBa/ii KOHLieHxpaL4iT bIahocho 
HOBMX flaHMX MoHxe-Kapjio (MK) fl/i5i npi/iMlxuBHUx MOflejiew ejieKxpojiixy +1: — 1 i +2:— 2. Bnpas Micxuxb 
xy nepeeary, u^o Hafljii/iLUKOBnti xiMinHi/iCi noxeHL4iaji MO>KHa oxpi/iMaxn npaiuio. He BflaKDHHCb flo xpyflOMicxKO- 
ro pospaxyHKy ri56ca-flK)reMa. Mw bw51bwjii/i, lu,o pesy/ibxaxi/i f/l Ha5iiw>KeHH!R fl,nff Koe4>ii4ieHxa cepeflHboT 
aKXMBHOcxi flo5pe y3rofl>KyioxbCiR 3 npoBefleHi/iMn fljia xaKoT caMoT MOflejii hmc/iobi/imh po3paxyHKaMW. KplM 
xoro, 6y/io npoxecxoBaHO xepMOflUHaMinHy y3rofl>KeHicxb f/l Ha5/iH>KeHH5i. Bi/ir/iaflae xaK, lu,o KoectDiL4ieH- 
xw cepeflHboT aKXHBHOcxi, nopaxoBahi 3a flonoMoroKD piBHjRHHiq ri65ca-flK)reMa, y3rofl>KyKixbciq 3 flaHi/iMi/i 
MoHxe-Kap/io xpoxw Kpame Hi>K BMpa3 raHceHa-Bi/iJie4)occa-Be/i/ioHi. fl/iiq L4ijiicHocxi po3paxyHKy npeflcxaB- 
neno xaKO>K fJl Hafljii/iLUKOBi BHyxpiiuHi eHepriT xa ocmoxumhI Koe43iL4ieHXM. U,\ pesy/ibxaxi/i nopiBHKDKDXbca 3 
po3paxyHKaMi/i Ha ochobI Ihlumx 3ara/ibHOBiflOMi/ix reopm, jiKi onucyioxb po3hhhm e/ieKxpo/iixiB, 30KpeMa 3i 

CepeflHbO-c4>epHHHHM Ha6/ll/l)KeHH!RM, MOflH4>iKaL4ieK> nixL4epa XeOpiT fle5a!R-rKDKKejl51, i rpaHllHHlIM 3aKOHOM 

fleGan-ftoKKejia. 

K/iKJHOBi c/iosa: npi/iMhi/iBHa MOfle/ib e/ieKTpojiiry, Koe4>ii4ieHT cepeflHbo'f aKTHBHOcri, hnepj}aHU,ioroBe 
Ha6jii/i>KeHHa, cepeflHbo-ccpepuHHe Ha6jii/i>KeHHa, MOfiejDOBaHHa MoHTe-Kapno, Merofl niTU,epa, Teopia 
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